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Topside Ionosphere 


ABSTRACT 


The existing uncertainties about the electron density profiles in the topside ionosphere, Le„ in the height regime from 
hmP2 to ->2000 km, requires the search for new data sources. Millions of io nograms had been recorded by the ISIS 
and Alouette satellites in the sixties and seventies, that never were analyzed in terms of electron density profiles. In 
recent years an effort started to digitize the analog recordings to prepare the io nog rams for computerized analysis. 
This paper shows how the digital ionograms are processed and the electron density profiles (from satellite orbit 
altitude, 1400 km for ISIS-2, down to the F peak) are calculated. 

The most difficult part of the task is the automatic scaling of the echo traces in the ISIS ionograms. 
Unlike the ionograms from modern ionosondes, the ISIS ionograms do not identify the wave polarization of the 
different echo traces, so physical logic must be applied to identify the ordinary (O) and extraordinary (X) traces, and 
this is not always successful. Characteristic resonance features seen in the topside ionograms occur at the gyro and 
plasma frequencies. An elaborate scheme was developed to identify these resonance frequencies in order to 
determine the local p la««n« gyrofreoaencies . This information helps in the identification of the O and X traces, 
and it provides the starting density of the electron density profile. 

The jflvtnton tf ttt wlrc Una iatt riwtrcn dauity nreDto «*** *•»« modified chebyshev 
polynomial fitting technique that is successfully applied in the ground-based Digisonde network. The automatic 
topside ionogram scaler with true height algorithm TOPIST is successfully scaling -70 % of the ionograms. An 
4 editing process* is available to manually scale the more difficult ionograms. 

The home page for the ISIS data project is at http;//ni>;>dc.ii>fc.M*>*.i;ov <suace/iMvisi>->taiu>.himl . It 
provides access to, as of January 2M1, 360.000 dig itized ISIS ionogram data and to related software. A search page 
lets users select data by location, time, and a host of other search criteria. The automated processing of the ISIS 
ionograms will begin later this year and the electron density profiles will be made available from the project home 
page. The ISIS data restoration efforts are supported through NASA's Applied Systems and Information Research 
Program. 
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Saving a Unique Data Set - The ISIS Story 


ISIS / launched in 1969 and /S/S 2 launched in 1971 were operated by NASA until 
1979. then by CRC (Canada) until 1984 and finally by RRL (Japan) until 1990. 
International participation - telemetry services, iongram reduction. 

Instruments: Topside Sounder, VLF, EPD, IMS, SPS, RPA, Photometers 

- Mayday - Mayday - (G. James, CRC) - The Candian Public Archive was getting 
ready to discard the more than 100 000 Alouette and ISIS telemetry tapes. 

ITM Data Evaluation Panel (D. Bilitza, Chair) selects the ISIS digitization project as 
its highest priority data restoration project. 

Some 13,800 Alouette and ISIS analog telemetry tapes were selected and shipped to 
GSFC. All but an additional - 4000 selected tapes have been destroyed. 

Processing begins; Analog tapes are read and digitized (R. Benson, PI. GSFC) and 
made online available fromlM ASA’s National Space Science Data Center (NSSDC). 

Topside lonogram Scaling and True Height Algorithm (TOPIST) software developed 
for automated scaling and inversion of ISIS ionograms. 

TOPIST processing of the more than 300,000 digitized ISIS-2 ionograms begins 

« Global mapping of the topside ionosphere over more than a solar cycle » 




Digitization of iSIS lonograms at GSFC 



System Used: The digitization system uses a 486-based Programmable Telemetry Processor (PTP), ao A/D converter 

board and a software device driver compatible with the OS2 operating system nsed by the PTP. Daring the A/D operation, 
ionogram frame sync pulses are identified, flagged aad the corresponding PB4 time code is determined. The PB4 time code is 
the output of a time-decoder card which receives the NASA 3*-b»t time code (added at the telemetry station) as the input The 
ionogram line sync pulses are also identified and flagged. These line-sync flags are used to identify the start of columns for 
data arrays used to generate digital lonograms. The flags are stripped off during the process of setting up these arrays. 


• RATEA/OLUME* Since processing started in 19% ISIS-2 iooograms have been di gitized 
at a rate of 5000 per month (7 GBytes) producing, as cf December 2000, a total of more 
than 300, 000 tanograms and a total data volume of 450 GBytes. 


Data Products: 

• Ionogram (3.75 km apparent-range (full) resolution), orbit parameter, OS2 binary 

The digitized sounder data consists of 8-bit receiver-amplitude values collected at a 40 kHz (0.025 ms) rate and the related echo 
tunes and frequencies. The frequency of each sounder pulse is determined from a 3rd degree interpolation between the times of the 
frequency markers as identified in the sounder video during the A/D operation. 

• Ionogram (15 km apparent-range (standard) resolution), orbit parameter, OS2 binary 

This data set is an average over four of the ‘full’ resolution data points. Typical size is 250-500 Kbyte, whereas the foil ionograms 
require typically 1-2 Mbytes. 

• Header files, oae per pass/tape, pass information and headers for all pass lonograms, ASCII 

• PCM, Pulse Code Modulation, one per pass, OS2 binary 

Consists of header, time. Automated Gain Control, frequency-marker information and data from all other instruments. 


Data, Soft ware, and Search Page 


Available Online from NSSDC: 

- All digitized ISIS ionograms (full and sta. res.), OS2 binary 

- All sta-res. ionograms in CDF format (created at NSSDC) 

- The header files, ASCII (PCM files are stored offline) 

- A Webinterface for searching for specific times, locations 
and other criteria (developed by R. Benson and team) 

- An interactive scaling anu inversion program (R. Benson) 

- The Worldmap program for computing orbit parameters 

- A selection of ionogram (about 1/10) is available for 
viewing on CDAWeb (plans for making all data available) 
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Total number of passes: 9,854 

Total number of ionograms: 302,809 
Total data volume (all files): 457 GB 

















Automatic Scaling of Traces 


Determine the possible trace points for each frequency line 
and get the the connection length to the nearest neighbors. 

Select the most likely points for inversion and identify the 
polarization (O/X) using an iteration technique. 

Select the best set of resonance and cutoff frequencies 

- Principle: Minimize the inversion error (error per point) for the largest 
possible set of trace points 

Fit generic function to trace 


! itlinti ueneric function u> trace 

R = R 0 - {[*o + DR\ exp[-a • {f-f s )\ - DR - exp [-6 • (f-f s )}} 
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FILE NAME: 7822206301 6KWA_A\A3 IS»S2TOPS_24S.OS2BIN 
SATELLITE: ISIS-2 ALTITUDE (km)-1 370.023 

LOCATION: GO Latttude/LongltucJe (Dog)— 10.787 / 165.764 
UTTIME: 06:30:16 Year/Day 1 978 / 222 Month/Day- 08/10 
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Electron density profile 


Auto-traces 


Ground echo 


Inversion of traces to e 

lec Iron 

density profiles 


The apparent range p 'for the frequencies f k assuming an electron density N(z): 


P ’(/*) = [m \N(z),f b fH(z), 0(2)] dz 

p group index of refractive media 
fH(z) gyrofrequencv 
0 dip angle 

Defining a reduced frequency: 

and new variable 

A for 0-trace 

f k ~ K fk ~fk * fH for X- trace 

1 fk ~fk * fH forZ-trace 

fk =M?r) 

2 f\-fN\z) 

fl-M 

Assume that Ne(h) can he represented with hehyshev polynomials: 


z ~ A m + g m % A ,r*(g), with g= \n(JNIfN m )l\n(JNJJN m ) 




Linear Equations for Coefficients 

The best coefficients in a least squares sense are 
found by minimizing the error 


A=iiP\- jam 

Linear equations for coefficients 

% A >Qij = t, p ’ k s jk 
Qtj = % s* % 
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scaling 


FILENAME: 731 21 184334BRZ_FUL_ISIS2TOPS 24S.OS2BIN 
SATELLITE: ISIS-2 ~ ALTITUDE (Kin)- 1300 

LOCATION: Geographic Latttuda/LongKucte (Deg)- -10.77 / 30.80 
UTTIME: 18:43:34 Y— r/Oay- 1073 / 121 Month/Oay- 05 / 01 
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TOPIST Performance 

• Test data: ~800 digital ISIS-2 ionograms 

• Current algorithm performance: 

- -70% of manually scalable ionograms are successfully autoscaled 

• Auto-scaling failures are caused by: 

- Severe spread 

- Data missing for a large frequency interval 

- Existence of unidentified traces (oblique?) 

- Incorrect identification of resonances 

• Inversion errors: 

- Ill-conditioned: No coefficients found 

- Non-monotonic profile found 

- N on-reasonable peak height found ( <200km) 

• Propagation modes: 

- The current TOPIST assumes that fH/f<l (for ISIS-2 sounding) and both O and X modes exist. The Z- 
trace is recalculated 


Future Plans 

• TOPIST testing and refining with a large 
number of ionograms from high and low 
latitudes. 

• Routine processing of the more than 
300,000 digitized ionograms into electron 
density profiles 





Treasure chest for ionospheric modeling 
- International Reference Ionosphere - 


•IRI topside profile is based on ~ 40,000 profiles from Alouette 1 

•Inversion of the now digitized ~ 300,000 iongrams will more 
than quintuple the current holdings of topside electron density 
profiles 

•Better ionospheric predictions for high solar activities. 

•Better understanding of the plasma processes in the high 
latitude and low latitude topside. 



